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ABSTRACT 

We estimate the expected yield of a non-targeted survey for galaxy clusters using the Sunyaev- 
Zel'dovich effect (SZE). Estimating survey yields requires a detailed model for both cluster properties 
and the survey strategy. We address this by making mock observations of galaxy clusters in cosmological 
hydrodynamical simulations. The mock observatory consists of an interferometric array of ten 2.5 m 
diameter telescopes, operating at a central frequency of 30 GHz with a bandwidth of 8 GHz. For a survey 
covering 1 deg 2 /month, we find that clusters with a mass above 2.5 x 10 14 H^qMq will be detected at 
any redshift, with the exact limit showing a modest redshift dependence. 

Using a Press-Schechter prescription for evolving the number densities of clusters with redshift, we 
determine that such a survey should find hundreds of galaxy clusters per year, many at high redshifts 
and relatively low mass - an important regime uniquely accessible to SZE surveys. Currently favored 
cosmological models predict ~ 25 clusters per square degree. 

Subject headings: cosmology: theory — galaxies: clusters: general — large-scale structure of universe 



1. INTRODUCTION 

The evolution of the cluster abundance is a sensitive 
probe of the mass density Q m (e.g 



1999, Bahcall and Fan 1998, Oukbir, Bartlett, and Blan- 



Viana and Liddk 



chard 1997 ). X-ray cluster surveys have started to con- 
strain Q m , but they are limited by their sample size 
and rapid decline in sensitivity with redshift, making the 
counts very sensitive to the selection function. While the 
selection function is presumably well-understood, it would 
be preferable to have a probe whose sensitivity does not 
fall off precipitously with redshift. We will show that a 
Sunyaev-Zel'dovich effect survey is ideal in this regard. 

Hot ionized cluster gas interacts with passing cosmic mi- 
crowave background (CMB) photons, distorting the CMB 
spectrum to create a decrement of CMB flux at lower fre- 
quencies and an excess at higher frequencies. This spec- 
tral distortion is independent of the redshift of the cluster, 
and only depends on the optical depth to Compton scat- 
tering and the temperature of the gas. This is the thermal 
Sunyaev-Zel'dovich effect (SZE) ( |Sunyaev and Zcl'dovicbl 
1972D . 

There have been several previous predictions of t he num- 
ber of clusters expected in SZE sur veys (e.g., Bartlett 
and Silk 1994, Barbosa et al. 1996), with most earlier 



work focusing on the total SZE flux. As survey yields de- 
pend sensitively on the observing strategy, we focu s here 
on yields fo r a proposed interferometric survey (see Mohi 
et al. 1999| ). 



form of curvature (open universe) or vacuum (cosmologi- 
cal constant) energy. Linear theory suggests that structure 
formation will slow considerably when the expansion dy- 
namics are no longer dominated by matter; this occurred 
around z <J 1.5 if current measurements of O m are correct. 
Therefore, the cluster abundance out to z ~ 2 should be 
a valuable probe of the matter density of the universe. 

A large collection of high-redshift, low-mass clusters also 
would provide an ideal sample for exploring evolution of 
the intra-cluster medium (ICM) and feedback from galaxy 
formation. The shallow potential wells of less massive clus- 
ters are more strongly affected by energy input from non- 
gravitational sources and are therefore the best place to 
search for the signatures of such processes. 

Estimating the cluster yield for an SZE survey requires 
that we know which properties of a cluster determine its 
likelihood of detection. We will show that the mass of a 
cluster is the single most important factor in determining 
if a given cluster can be detected by an interferometric 
survey. In this case, the calculation of the expected yield 
separates into two distinct exercises: finding the minimum 
observable mass as a function of redshift and calculating 
the number density of clusters above a given mass thresh- 
old as a function of redshift. 

We determine the minimum observable mass as a func- 
tion of redshift by making synthetic observations of N- 
body+gas hydrodynamical simulations (§0). The number 
density of cl usters is calculated using the Press-Schechter 
prescription (Press and Schcchtcr 1974| ; §g). Section | con- 
tains estimates of the expected survey yield and the results 
are discussed in Section [|. 



A large catalog of high-redshift clusters that extends 
to low masses would be an extremely useful resource for 
several reasons. Observations suggest that the universe 
may no longer be matter-dominated, with a large frac- 
tion (~ 70%) of its present energy density either in the 
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2. CLUSTER DETECTABILITY 

The SZE decrement along a given line-of-sight is inde- 
pendent of cosmology. It is simply proportional to the inte- 
grated thermal pressure along the line-of-sight and there- 
fore only depends on the properties of the cluster. The 
decrement can be written 



=g[x) / dln e {l) —<J T , 

J-CMB J m e c z 



where n e is the electron number density, T e is the elec- 
tron temperature, m e is the electron rest mass, <jt is the 
Thomson cross section, g(x) = x(e x + l)/(e x — 1) — 4 with 
x = hv /ksTcMB and the integral is along the entire line- 
of-sight; by assumption, the only significant contribution 
to the integral comes from the cluster atmosphere. In the 
Rayleigh- Jeans limit [y <C 200 GHz) the dimcnsionless 
frequency factor g(x) = — 2. 

The specific intensity S v in the Rayleigh- Jeans regime is 
S u — 2kB&Tv 2 1 'c 2 d£l, where d£l is the effective solid an- 
gle of the observations. Thus the total SZE flux decrement 
Stot for a galaxy cluster can be written 

2k 2 B v 2 g(x)a T TcMB /rr , M 2Q ofiCM 

Otot(Z) - 4 , , x 2 {J-e/n > \ Z > 



m e c 4 dA(z) 



where dA{z) is the angular diameter distance, (T e ) n is the 
electron density weighted mean temperature in the clus- 
ter, fi e is the mean molecular weight per electron, m p is 
the proton mass, ficM is the ratio of total gas mass to 
binding mass, and M200 is a measure of the cluster virial 
mass, defined as the mass within r2oo, the radius where 
the mean interior density is 200 times the critical density. 
Note that we explicitly ignore contributions to the SZE 
flux coming from outside the virial region. 

Equation 2 indicates that the total SZE flux for a clus- 
ter is directly proportional to the cluster virial mass; the 
only dependence on cluster structure is through the den- 
sity weighted mean temperature (T e ) . This is unique to 
a survey with a beam large enough that the cluster SZE is 
not resolved. However, for maximum brightness sensitiv- 
ity to the SZE effect, the beam used for a survey should 
be well matched to the typical angular scale subtended by 
clusters. Such a survey will partially resolve most of the 
clusters and will therefore be somewhat sensitive to the 
internal cluster structure. We show below that for a pop- 
ulation of clusters with similar ICM mass fraction ficM 
and similar temperature structure, the detection threshold 
for an interferometric SZE survey is also effectively a virial 
mass limit. 

2.1. Mock Interferometric SZE Observations 

We determine the detection threshold of an SZE survey 
by analyzing mock observations of numerical cluster sim- 
ulations. These mock observations are appropriate for a 
proposed ten-element interferometer, composed of 2.5 m 
diameter telescopes outfitted with receivers operating at a 
central frequency of 30 GHz and a bandwidth of 8 GHz. 
We assume a correlator efficiency of 0.88 and an aperture 
efficiency of 0.77. We take the system temperature below 
the atmosphere to be 21 K; this includes contributions 
from spillover past the primary. We assume an atmo- 
spheric column with opacity r = 0.045 at zenith, which 



is a conservative estimate appropriate for a low altitude, 
moderately dry site such as the Owens Valley Radio Ob- 
servatory in summer. The integration time on each cluster 
is 42 hr, composed of six 7 hr tracks on the cluster. This 
exposure on each piece of the sky would allow us to survey 
~1 deg 2 per month. 

Interferometers measure the visibility, V(u, v) which is 
the Fourier transform of the sky brightness distribution 
multiplied by the primary beam of the telescopes. We cre- 
ate mock observations at a redshift z by imaging hydro- 
dynamical cluster simulations (described in detail below) 
at that redshift; we place those clusters at the appropriate 
angular diameter distance, multiply the resulting SZE im- 
age with the primary beam of the 2.5 m dishes (modeled 
as a Gaussian with FWHM = 14.7' at 30 GHz), and then 
Fourier transform to produce visibilities V(u,v). We then 
sample these visibilities at the same locations in u—v space 
which appear in our simulated 7 hr array track and add 
the appropriate noise. 

2.2. Hydrodynamical Cluster Simulations 

The effects of ongoing cluster merging on the ICM den- 
sity and temperature structure can be calculated self- 
consistently in hydrodynamical simulations. Therefore, 
these simulations provide a way of producing test clus- 
ters whose complexity approaches that of observed clus- 
ters. In this work we use an ensemble of 36 hydrody- 
namical cluster simulations carried out within three dif- 
ferent cold dark matter (CDM) dominated cosmologies (1) 
SCDM (n m = 1, a s = 0.6, h 50 = 1, T = 0.5), (2) OCDM 
(fi m = 0.3, cr 8 = 1.0, h 50 = 1.6, T = 0.24), and (3) LCDM 
(fi m = 0.3, Q A = 0.7, a 8 = 1.0, h 50 = 1.6, T = 0.24). 
Here erg is the power spectrum normalization on 8/i _1 Mpc 
scales; initial conditions are Gaussian random fields con- 
sis tent with a CDM transfer function with the specified 
r (Davis et al. 1985). These simulations have been used 
previousl y for studies of the X-ray emission from galaxy 
clusters ([Vlohr and Evrard 1997 , Mohr, Mathiescn, and| 
Evrard 1999^ 

Within each cosmological model, we use two 128 3 N- 
body only simulations of cubic regions with scale ~ 
400 Mpc to determine sites of cluster formation. Within 
these initial runs the virial regions of clusters with Coma- 
like masses of 10 15 M Q contain ~10 3 particles. 

Using the N-body results for each model, we choose 
clusters for additional study, resimulating them at higher 
resolution with gas dynamics and gravity, as described 
below. The size of the resimulated region is set by the 
turnaround radius for the enclosed cluster mass at the 
present epoch. The large wavelength modes of the ini- 
tial density field are sampled from the initial conditions of 
the large scale N-body simulations, and power on smaller 
scales is sampled from the appropriate CD M power spec - 
trum. The simulation scheme is P3MSPH flEvrard 1988| ), 
the ICM mass fraction is fixed at ficM = 0.2, and radia- 
tive cooling and heat conduction are ignored. 

The high resolution, hydrodynamical simulations of in- 
dividual clusters require two steps: (1) an initial, 32 3 , 
purely N-body simulation to identify which portions of 
the initial density field lie within the cluster virial region 
at the present epoch, and (2) a final 64 3 , three species, 
hydrodynamical simulation. In the final simulation, the 
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portion of the initial density field which ends up within 
the cluster virial region by the present epoch is represented 
using dark matter and gas particles of equal number (with 
mass ratio 4:1), while the portions of the initial density 
field that do not end up within the cluster virial region 
by the present epoch are represented using a third, colli- 
sionless, high mass, species. The high mass species is 8 
times more massive than the dark matter particles in the 
central, high resolution region. This approach allows us 
to include the tidal effects of the surrounding large scale 
structure and the gas dynamics of the cluster virial region 
with simulations that take only a few days of CPU time 
on a low end UltraSparc. 

The scale of the simulated region surrounding each clus- 

1 /3 

ter is in the range 50-100 Mpc, and varies as M h ' alol 
where M^aio is approximately the mass enclosed within 
the present epoch turnaround radius. Thus, the 36 sim- 
ulated clusters in our final sample have similar fractional 
mass resolution; the spatial resolution varies from 125- 
250 kpc. We will find that the clusters of greatest interest 
to us are those with mass ;> 2 x 10 14 M Q , which have at 
least several thousand gas particles, even in the lowest res- 
olution simulations (i.e., the highest mass at z = 0). The 
masses of the final cluster sample vary by an order of mag- 
nitude. Following procedures described in Evrard (1990), 
we create SZE decrement images along three orthogonal 
lines of sight for each cluster. Each image is 128 2 and 
spans a distance of 6.7/i^q 1 Mpc in the cluster rest frame. 

A strength of using numerical cluster simulations is 
that structural evolution consistent with the cosmologi- 
cal model is accounted for naturally by simply examining 
higher redshift outputs of the simulations. 
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Fig. 1. — Detection significance, Ax 2 relative to null model, as 
a function of cluster mass for an interferometric SZE survey. All 
clusters were imaged at a distance corresponding to z=0.5 in an 
f2m = 0.3 and Q\ = 0.7 universe. The dashed horizontal line corre- 
sponds to a 5cr detection for a two-parameter fit. 



2.3. Determining the Survey Detection Threshold 



We attempt to detect the clusters in the mock obser- 
vations by fitting th e data to the Fourier transforms o f 
spherical (3 models ( Cavaliere and Fusco-Femiano 1978 ), 
with (3 — 4/3. The results are not sensitive to the value 
of j3, and in fact a simple Gaussian works well. We choose 
this value of (3 because its transform is a simple exponen- 
tial and because the best fits to the SZE decrement images 
from the simulations yield a mean value (/?) = 1.1. We de- 
termine the best fit core radius 9 C and central decrement 
AT by minimizing \ 2 ; the Ax 2 difference between the 
best fit (3 model and the null model (no cluster present) 
provides a measure of the significance of the cluster detec- 
tion. We set our threshold Ax 2 > 28.7 (5u for two degrees 
of freedom) . 

Figure [j] shows the detection significance for 108 mock 
observations of the simulated clusters output at redshift 
z = 0.5. Observations of clusters in all three cosmologies 
appear on the same plot, and for this example we have im- 
aged all the clusters at the same angular diameter distance 
and accounted for differences in H a . There is a striking 
correlation between the detection significance and cluster 
virial mass, indicating that even when complex cluster dy- 
namics are considered, the survey detection threshold can 
still be effectively described as a mass threshold. The scat- 
ter about this correlation is a reflection of the variation in 
cluster structure due to different merger histories and pro- 
jection effects. We determine the mass threshold by ex- 
amining the Ax 2 -M2oo relation and determining the virial 
mass at which Ax 2 = 28.7. We calculate the RMS scatter 
about the relation at A^ 2 > 28.7, and use that scatter 
as an estimate of the width of the detection threshold in 
mass; the number of mock observations used to character- 
ize the RMS varies between 21 in our highest redshift bin 
to over 100 in our low redshift bins. We model the scatter 
as a Gaussian distribution in x 2 a t each mass. In this way, 
the fraction of detected clusters at each mass is expressed 
as an integral over a Gaussian. 

Note that the small scatter about the Ax 2 -M 2 oo relation 
indicates that mass is indeed the primary factor in deter- 
mining cluster detectability. Moreover, the fact that all 
three cosmologies produce consistent relations once differ- 
ences in H a are accounted for indicates the relative insen- 
sitivity of the detection threshold to differences in cluster 
structure. 

By repeating this exercise at several redshifts, we deter- 
mine the survey mass threshold as a function of redshift, 
shown in Figure ^. Note that the thresholds differ for each 
cosmology because of differences in the angular diameter 
distance- redshift relation oIa{z). The error bars on the 
mass limit points indicate the RMS scatter in mass about 
the Ax 2 - M 2 oa relations (see Fig. |l|). 

The hydrodynamical simulations cannot be used to ex- 
tract mass thresholds beyond redshift z ~ 2.3, because no 
clusters in our ensemble are massive enough to lie above 
the detection threshold. For higher redshift and to aid in 
interpolating the mass threshold at arbitrary redshift, we 
use mock observations of spherical (3 models normalized to 
agree with the simulations (curves in Figure |^) . We evolve 
the /3-models in redshift using the spherical collapse model 
(Lahav et al. 1991); this evolution is self-similar, account- 
ing for the difference in cluster structure due to evolution 
of the mean cosmological density. We use these smooth 
curves essentially as fitting functions. 
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The insensitivity to redshift of the minimum detectable 
mass in an interferometric survey follows from a balanc- 
ing of several effects. The increasing angular diameter 
distance cLa with redshift decreases the total cluster SZE 
flux as d^ 2 (see equation 2), tending to increase our lim- 
iting mass. However, this effect is largely offset by cluster 
evolution. At higher redshifts clusters are denser, and, at 
constant virial mass J Vhno, have higher virial temperatures 



where p is the mean comoving background density, a(M, z) 
is the variance of the fluctuation spectrum filtered on mass 
scale M, and 6 C is the effective linear overdensity of a per- 
turbation which has collapsed and virialized. In principle, 
5 C has a modest dependence on the cosmological density 
parameter; the spherical collapse model predicts a varia- 
tion of only ~ 5% for fi m ~ 0.1 — 1. In this work , we 



T. Beth of these effects increase the total SZE flux. In ad- 



assume a constant threshold S c = 1.69 for simplicity (Pee- 



bles 1980) 



dition, at higher redshift a cluster has a smaller apparent 
size, which enhances the cluster visibility V(u, v) (Fourier 
transform of SZE decrement distribution; see §2.1) at the 
baselines where we make our measurements. Taken to- 
gether, these effects largely explain the behavior of the 
limiting mass with redshift. 



Following Viana and Liddle (1999), we take the variance 
in spheres of radius R to be 



a(R,z) = ag(z)( 



R 



Sh-^Mpc 



where 



o 



1 r 

0.9 I- 
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7 (i?) = (0.3r + 0.2) 



2.92 + log! 



-7(R) 



R 



&h~ l Mpc 



(5) 



(6) 
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The comoving radius R is determined as the radius which 
contains mass M at the current epoch, while T is the usual 
CDM shape parameter, taken to be 0.25 for this study 
(Peacock and Dodds 1994, Dodelson and Gaztanaga 1999) 
unless stated otherwise; we show that the results are in- 
sensitive to the exact choice of T. 

We examine the cluster abundance in three cosmolog- 
ical models: oCDM (O m = 0.3, n A = 0, h = 0.65, T = 
0.25, ct 8 = 1.0), ACDM (fi m = 0.3, H A = 0.7, h = 
0.65, T = 0.25, cr 8 = 1.0), and rCDM (fi m = l,fi A = 
0, h = 0.5, r = 0.25,(78 = 0.56), with the last model sim- 
ply a CDM model with the transfer function modified to 
agree with observations of galaxy clustering. Note that 
these models differ slightly from the cosmologies assumed 
for the simulations. 

All models are chosen to have a global ICM fraction 
fiCM = 0.2, in rough agr eement with observed ICM 
mass fractions of clusters (David, Jones, and Forman 



betwee i cosmologies are simply showing the cosmological depen- 
dence ot tne angular diameter distance— rcdshilt relation, f rom top 
to bottom, results are shown for ACDM (solid triangles and dot- 
ted line), oCDM (open triangles and dashed line) and tCDM (solid 
squares and solid line). The uncertainties in the mass threshold are 
the same for all cosmologies but, for clarity, are only shown for the 
tCDM points. 

THE CLUSTER ABUNDANCE AND ITS REDSHIFT 



19951, IWhitc and Fabian 19951 , Grcgo et at 200q , Mohr 



EVOLUTION 



Mathicscn, and Evrard 1999). As can be seen from cqua- 
tion 2, the mass limits will depend on the ICM mass frac- 
tion and therefore the expected yields will also be sensitive 

to fiCM- 

We show that the expected survey yield is very sensi- 
tive to as, also calculating the expecte d yields for oCDM 
with a lower value of erg = 0.85 (e.g., Viana and Liddle 
1999). The constraints on erg will improve dramatically 



Given a minimum observable mass, one can calculate 
the number of observed clusters as: 



N(M > M thres ) 

cfdnjdzj^ thr 



,(*) 



dM 



dn(M,z) d A (z) 2 (l+z) 2 



dM 



(3) 



where eifi is the solid angle and n{M, z) is the comoving 
number density. To calculate the comoving number den- 
sity of clusters, we used th e Press-Schechter prescription 
(Press and Schechter 1974). 

The comoving number density of bound objects between 
masses M and M + dM at redshift z is given by 



in the near future, as new X-ray telescopes are expected 
to provide a much better determination of the local abun- 
dance, so we do not expect uncertainties in ers to affect 
interpretation of survey results. 

In a critical density universe (fi m = 1,Oa = 0), ers °c 
(l + z) -1 . Following Carroll, Press, and Turner (1992), we 
express growth in alternate cosmologies through a growth 
suppression factor which can be approximated as 



1 



fi m /2) (l + r> A /70 



dn(M, z) 
dM 



2 p da{M,z) 



tt M dM a 2 (M,z) 



exp 



2ct 2 (M, z) 



In this notation, we can now express the normalization of 
the power spectrum as 



(4) 



08 0*) = 



,(0) g(Q m (z),n A (z)) 



1 + z g(Q m (0),Q A (0)) 



(8) 
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4. EXPECTED SURVEY YIELD 

Figure || shows both the differential counts as a func- 
tion of redshift and integrated number of clusters for our 
three cosmologies. Cluster physics is especially uncertain 
at high redshift, so we have chosen to cut off all integrals 
at a redshift of z = 4. 
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Fig. 3. — Differential counts (top panel) and integrated counts 
(bottom) as a function of redshift per square degree for two differ- 
ent methods of dealing with the limiting mass. The solid curve for 
each cosmology assumes a step function at the best fit limit, while 
the dotted curve corresponds to a Gaussian distribution of \ 2 a t 
each mass. 

These results are quite exciting, because they indicate 
that an SZE survey will yield a large cluster catalog ex- 
tending to high redshift. If the currently favored ACDM 
model is correct, then we expect to detect 300 clusters in 
a one year survey. Figure |^ also shows that the width of 
the mass limit has little effect on the cluster yield. The 
yield for each cosmology is shown with two curves; the 
solid curve corresponds to a step function mass limit where 
A\ 2 = 28.7, and the dashed curve corresponds to a mass 
limit modeled as a Gaussian probability distribution in 
Ax 2 at each mass, with variance set by the scatter around 
the Ax 2 — M200 relation seen in the mock observations. 
In particular, the largest differences arise mainly in the 
high-redshift region, where there are few clusters above our 
mass threshold. While encouraging, this threshold uncer- 
tainty requires more attention before future SZE surveys 
can be correctly interpreted. 
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Fig. 4. — Differential (top panel) and integrated (bottom) ex- 
pected number of cluster in an oCDM model for two values of the 
normalization of the power spectrum, erg = 0.85 (dotted) and 1.0 
(solid). 

The differences between cosmologies at z > 1 are very 
promising. Differences arise because of the different rates 
of structure growth, the dependence of the mass threshold 
on the angular diameter distance, and the differences in 
the volume element of the survey. Open models have both 
a smaller angular diameter distance and more structure 
at high redshift, both leading to more observed clusters. 
The effect of ACDM probing a larger volume is a relatively 
small effect. 

The expected cluster yield depends on cosmology 
through the growth rate of perturbations (fixed mainly by 
the density parameter f2 m ), the volume per unit solid an- 
gle and redshift, and the shape and normalization of the 
initial power spectrum. Varying T, the CDM shape pa- 
rameter, within the 95% confidence interval, 0.2 < V < 0. 3 
(Peacock and Dodds 1994, Dodclson and Gaztanaga 1999), 
leads to increases (T = 0.3) or decreases (T = 0.2) in the 
integrated counts of 6% for ACDM, 8% for oCDM and 
17% for tCDM. The effect is small, indicating that SZE 
survey counts alone will not strongly constrain the shape 
of the power spectrum. However, cluster samples con- 
structed through SZE surveys can be used to study the 
shape of the power spectrum through other means, such 
as the two-point correlation function. 

On the other hand, the normalization of the power spec- 
trum, parameterized by as, is quite important for pre- 
dicted SZE counts, as demonstrated in Figure ||. For clar- 
ity, we have only shown the effect for oCDM. The uncer- 
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tainty in eg results in a large uncertainty in the cluster 
yield. Upcoming X-ray observations of nearby clusters are 
expected to constrain as to much higher accuracy, and the 
proposed SZE survey will provide an independent mea- 
surement. Thus, while the erg uncertainty is significant for 
predicting yields, it will not be a serious impediment to 
interpreting results from an SZE survey. It is important 
to no te that, even for low values of erg, we expect to find 



a significant number of high-rcdshift clusters. 



crease possible CMB contamination significantly. Indeed, 
larger telescopes may even be slightly more efficient for 
detecting high-redshift low-mass clusters, as the smaller 
primary beam is better matched to the compact nature of 
these objects. 

We use an ICM mass fraction of 20% in these simu- 
lations, which is inconsiste nt with the observational con- 
straints, ficM = 0.21/i^n 1,5 (Mohr, Mathiesen, and Evrard 
Grego et al. 2000), if H is significantly higher 
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5. DISCUSSION 

In a low-density universe, half of the clusters in an SZE 
survey will have z J> 1. While other surveys may yield 
comparable numbers of clusters at high redshift by sur- 
veying a larger fraction of the sky, the SZE catalog will be 
unique in having similar sensitivity at all redshifts; SZE 
surveys are effectively limited only by the abundance of 
clusters above the lowest observable mass. PLANCK is 
expected to find ~ 10 4 clusters, but the effective limiting 
mass is f airly high, resulting in relatively few high-rcdshift 



cluste rs (De Luca, Desert, and Puget 1995, da Silva et al 
1999D 



An interferometric survey with a synthesized beam that 
partially resolves clusters has several advantages. Point 
sources are easily identified and removed from the data; 
thus, it is unlikely that point sources will systematically 
affect the magnitude of the observed cluster decrement. 
In addition, the more massive clusters in our sample can 
be imaged with high S/N at the same time as they are 
detected. 



than 50 kms _1 Mpc _1 . However, simple virial arguments 
(T cx M 2 / 3 ) and Eqn. ^ indicate that the limiting mass 
scales as fjQM- For H° = 80kms _1 Mpc _1 , we expect 
fiCM = 0.10. While a gas fraction of only 10% would se- 
riously reduce the expected number of clusters (by a fac- 
tor of 2.5-3) the resulting catalog would still be large and 
unique. 

An uncertainty that we have not discussed is the effect 
of galaxy formation or preheating on the ICM structure 
and SZE detectability. While this is being addressed with 
numerical simulations currently in progress, we do not be- 
lieve that our results should depend significantly on gas 
evolution. ICM evolution mainly affects the core regions 
of clusters ( Ponman, Cannon, and Navarro 1999| ), on an- 
gular scales smaller than those for which the brightness 
sensitivity of the interferometer is optimized, 2' to 7'. The 
total SZE flux from the unresolved core depends only on 
the temperature and the number of electrons at that tem- 
perature (see §H) . Even fairly extreme models of gas evo- 
lution lead to only small changes in the expected counts 
for a survey sensitiv e primarily to total SZE flux ( [Holdci 



Ex pensive follow-up will be required to make best use and Carlstrom 1999| ) and we expect this to be true for any 



of an SZE survey. Extensive optical observations will be 
required to obtain redshifts for the detected clusters. This 
may be difficult for the highest-redshift objects, but it is 
not an insurmountable problem. Follow-up with X-ray 
telescopes would be helpful as well, but the low-mass, 
high-redshift objects are expected to be undetectable in 
the X-ray band, even with 10 5 s XMM exposures. How- 
ever, redshifts, ICM temperatures and X-ray images for a 
portion of the sample would enable direct SZE+X-ray dis- 
tances, such as those being meas ured with the currently 
available data (Reese et al. 1999 and references therein). 
These distances constrain the angular diameter distance 
relation, which provides an independent measurement of 
the cosmological matter density; these results would be 
complementary to those available from consideration of 
the cluster counts alone. 

Deep exposures such as the ones considered here could 
be contaminated by primary anisotropies in the CMB. 
A minimum separation of 2.5m corresponds to multipole 
I = 1571 at 30 GHz. At this scale, the CMB anisotropy 
levels could well be larger than 10fiK, which would be 
above the noise levels that we have assumed but well be- 
low the detection threshold. A shift to higher t may be 
required to avoid these effects, which can be achieved by 
longer baselines (which could allow larger telescopes) or a 
higher observing frequency. We have found that a shift to 
I ~ 3000 can be easily accommodated with only a small 
loss in cluster detection efficiency. The most numerous 
clusters will be the low-mass clusters near the detection 
threshold, which will be the most compact. Because of 
this, moving to higher multipoles would not severely af- 
fect our sensitivity to these objects, while this would de- 



survey that does not resolve the core regions of clusters. 

While the expected cluster yield should be fairly insensi- 
tive to ICM evolution, the properties of the cluster sample 
and their evolution with redshift will shed considerable 
light on the question of evolution of the ICM. At the same 
time, such a survey will provide determinations of key cos- 
mological parameters that will be entirely independent of 
all other determinations. 
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